say asap 


=< AFAPL~TR-78-6 al 
Part III 
t- t. eS 
> aS - 
= p° 
xO 
ROTOR-BEARING DYNAMICS TECHNOLOGY DESIGN GUIDE 
Part III 
Tapered Roller Bearings 
A. B. Jones 
J. M. McGrew, Jr. 
Shaker Research Corporation 
Northway 10 Executive Park 
Ballston Lake, New York 12019 
February 1979 
> 
Chun 
) 
jo Interim Report for Period April 1977 ~ December 1978 
Lu ‘ 
= 
ue Approved for public release; distribution unlimited. 
\ c.? 
Cc 
ce 


AIR FORCE AERO PROPULSION LABORATORY 


AIR FORCE WRIGHT AERONAUTICAL LABORATORIES 
AIR FORCE SYSTEMS COMMAND 


WRIGHT-PATTERSON AIR FORCE BASE, OHIO 45433 


REPRODUCED FROM 


BEST AVAILABLE COPY 79 | 6 | 0 6 2 4 


NOTICE 


When Government drawings, specifications, or other data are used for any 
purpose other than in connection with a definitely related Government pro- 
curement operation, the United States Government thereby incurs no responsi- 
bility nor any obligation whatsoever; and the fact that the Government may 
have formulated, furnished, or in any way supplied the said drawings, speci- 
fications, or other data, is not to be regarded by implication or otherwise 
as in any manner licensing the holder or any other person or corporation, or 
conveying any rights or permission to manufacture, use, or sell any patented 
invention that may in any way be related thereto. 


This report has been reviewed by the Information Office (OI) and is releas- 
able to the National Technical Information Service (NTIS). At NTIS, it will 
be available to the general public, including foreign nations. 


This technical report has been reviewed and is approved for publication. 
} 


B. 
JOHN B. SCHRAND 
Project Engineer 


FOR THE COMMANDER 


7 et 
BLACKWELL C, DUNNAM 
Chief, Fuels and Lubrication Division 


If your address has changed, if you wish to be removed from our mailing list, 
or if the addressee is no longer employed by your organization please notify 
AFAPL/SFL, WPAFB, OH 45433 to help us maintain a current mailing list. 


Copies of this report should not be returned unless return is required by 


security considerations, contractual obligations, or notice on a specific 
document. 


AIR FORCE/56780/23 August 1979 — 140 


SECURITY CLASSIFICATION OF THIS PAGE (When Date Entered) 
¢ 


“(0.7 {REPORT DOCUMENTATION PAGE a 
2 12. GOVT ACCESSION NOW 3. RECIPIENT'S (o> NUMBER 
wupmere-oerece ne T=3) | 


Ta TITLE (and Subtitte) BT YRE OT REROR A ; 

save Le eee ts ws crac re Technical ail ate ye 
_ROTOR~BEARING DYNAMICS. TECHNOLOGY DESI ; 

tee aia GN (GUIDE “iP A rid 1977 = — November 1978 
Tapered Roller Bearings « ‘ e 


7. AU THOR). 
ee 


ee Psat on OO eau ne OZ 


A. B. fiouva f 
J. M. /McCrew, Jr. 


; ¥33615-76-C-2038 | 
9. PERFORMING ORGANIZATION NAME AND AOORESS 10. PROGRAM ELEMENT, PROJECT, TASK 
AREA & WORK UNIT NUMBERS 

Shaker Research Corporation Aor f Or 

Northway 10 Executive Park & es 

Ballston Lake, N.Y. 12019 3048-6 a8 
18. CONTROLLING OF FICE NAME ANO AODRESS % 7 
Air Force Aero Propulsion Laboratory/SFL 
Air Force Systems Command 


Wright-Patterson AFB, 0 
. MONITORING AGENCY NAME & AOCORESS/(if difterent from Controlling Office) 18 SECURITY CLASS. (of thte report) 


CY “ Unclassified 
i Se. CECL ASSIFICATION/ COWNGRADING 
SCHEDULE 


Approved for public release; distribution unlimited. 


e 


80 
16. DISTRIBUTION STATEMENT (of thie Report) 


17. DISTRIBUTION STATEMENT (of the ebatract entered in Block 20, !f different trem Report) 


18. SUPPLEMENTARY NOTES 


19. KEY WORDS (Continue on reveree side if necessary and identity by block mamber) 


Tapered Roller Bearings Roller Bearing Stiffness 
Tapered Roller Bearing Stiffness Turbine Bearings 
Roller Bearings Rotordynamics 


20. A@STRACT (Continue on reverse side if neceseary and identity by bieck mamber) 

This report is an update of the original Part IV of the Rotor-Bearing 
Dynamics Design Technology Series, AFAPL-TR-65-45 (Parts I through X). 

A computer program is given for preparation of tapered roller bearing 
stiffness data input for rotordynamic response programs. The complete 
stiffness matrix is calculated including centrifugal effects. Consider- 
ations such as elastohydrodynamic and cage effects are not included since 
chey have eaters eb ence on the catmaiat ion of saa aps roller bear ios 


DO , rom 1473s EDITION OF | wov 68 18 OBSOLETE \ “3 
S/N 0102-014- 6601 | : 
SECURITY CLASSIFICATION OF THIS PAGE (When Date Batwod) (4 


oO #. fy 
SOL # ee. y cia 


FOREWORD 


This report was prepared by Shaker Research Corporation under USAF 

- Contract No. AF33615-76-C-2038. The contract was initiated under Project 
3048, "Fuels, Lubrication, and Fire Protection," Task 304806, "Aerospace 
Lubrication," Work Unit 30480685, "“Rotor-Bearing Dynamics Design." 


The work reported herein was performed during the period 15 April 
1977 to 15 November 1978, under the direction of Jciin B. Schrand (AFAPL/ 
SFL) and Dr. James F. Dill (AFAPL/SFL), Project Engineers. The report 
was released by the authors in December 1978. 


iii 


Section I 


Section II 


Section III APPLICATION OF COMPUTER PROGRAM 


TABLE OF CONTENTS 


INTRODUCTION 


ANALYSIS 


wee ew mee ww weeeeweowecocewe ses ww eet eer eccwoo cece 


2.1 General Bearing Model and Coordinate System 


2.2 General Bearing Support Characteristics 


eewooewre 


2.3 Tapered Roller Bearing Characterization 


er 


2.4 Tapered Roller Bearing Under Combined Loading _ 


3.1 Sample Test Case 


wwe ee mem wee ewe wee meme ne ereowweee 


3.2 Input Format 


www ee we wm mmm ewe we we wee ee wee weer oeecenee 


3.3 Output Format 


ewww wwe nw we we ewe ww we wm ee wwe emer 


APPENDIX COMPUTER PROGRAM FOR CALCULATING STIFFNESS MATRIX 
OF TAPERED ROLLER BEARING. == eo 
REP ERENCES Nahin ie Sue aia ee ana net Gets ey seon Owes 
Vv 
PRECEDING PAGE BLANK 


an ig 


noe 


oat ie es 


80 


Figure 
Figure 
Figure 
Figure 
Figure 
Figure 
Figure 
Figure 
Figure 
Figure 
Figure 
Figure 
Figure 
Figure 
Figure 


Figure 


LIST OF ILLUS TRATIONS 


Bearing Stiffness Model 


we ee eee oe a 2 eee cesceeeocee en ewee 


Bearing Lecation Coordinate System 


Linearization of Tapered Roller Bearing Stiffness 


Tapered Roller Bearing 


Bearing Coordinate System 


Tapered Roller Bearing Index, q 


Dimensions of Crown Drop 


Forces and Moments on Roller 


Output Data for Load Condition #1 


ew mew mowe eee e eo ewes. 


Output Data for Load Condition #2 


vi 


45 


49 


(8) 1 ineal 
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NOMENCLATURE 


Description 


Semi-width of contact ellipse at x 
Corner break at roller small end 


Corner break at roller big end 


Damping component, change of force in i direction 


due to velocity in j direction; i = x, y, 23 
j =X, y, 2. 


Damping matrix 


(By) 1ineal 0 0 

0 0 
0 0 

(By) angular 
0 0 


Damping matrix due to lateral velocities 


B " B 
y yy ineal 


Damping matrix due to angular velocities 


B B 
XX xy 
Bs B 
y yy ngular 
l fori=l 
A constant, Cc, = { 
-l for i = 2 


Roller diameter at midpoint of effective length 
of roller 


Roller diameter at x 
Pitch diameter at midpoint of effective length 


Modulus of elasticity for roller 


vil 
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Units 
in. 
in. 
in. 

lb-sec 


in 


lb-sec 
in 


ib=sec 


in 


in-lb-sec 


radian 


in. 


tha/in” 


I 


(KE) ineal 


Modulus of elasticity for race body ibs /in* 


Pitch diameter at x in. 
Roller centrifugal force lbs. 
External applied force, i= x, y, z lbs. 
Reaction force, positive in direction opposite lbs. 


to displacements, i x, y, z 


F 
Force Matrix = ae "Ibs. 
F 
Zz 
Distance along roller cone element from extreme in. 


end of effective length to point where crown drop 
is measured 


Roller crown radius minus the rise of the arc at in. 
midpoint of effective length 


Moment of inertia about roller center of gravity ibs-in® 
Roller-race stiffness Ibs/in. 
Stiffness component, change of force in i direction lbs/in. 
due to displacement in j direction. 

i =x, y, 2 j™ Xx, y, z 


Stiffness matrix 


(KY 1ineal 9 _ 


0 0 
0 0 
(Ky) angular 
0 0 
Stiffness matrix due to lateral displacements lbs/in. 
K K 
xx xy 
K ‘ K 
y Y¥4} ineal 
Stiffness matrix due to angular rotations in-1b 
K K tad 
xX xy 
K K 


yx yy angular 


viii 


Perpendicular distance from the line of action 
of flange reaction, P3, to roller centerline at 
midpoint of effective length 

Effective length of roller load carrying surface 


Length of flat portion of roller measured along 
roller cone element 


Total length of roller measured parallel to roller 
axis between sharp intersections of end faces with 
roller cone elements 

Mass of roller 

. Jller gyroscopic moment 

External applied moment, i = x, y, z 
Reaction moment, i = x, y, z 

Outer race/roller contact moment 

Inner race/roller contact moment 

Number of rollers 

Outer ring rotational speed 

Inner ring rotational speed 

Contact unit loading 

Current estimate of contact unit loading 
Diametral clearance 

Outer contact load on qth roller 

Inner contact load on qth roller 

Flange reaction on qth roller 

Roller position index 

Rolier crown radius 

Roller big-end spherical radius 


Radial distance from roller center line to 
flange reaction 
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in. 


in. 


in. 


in. 


lbs. 
lbs-in. 
lbs-in. 
lbs-in. 
lbs-in. 


lbs-in. 


rad/sec. 
rad/sec. 


lbs/in. 


lbs/in. 


Bearing coordinate system 


tatic component of displacement 
Dynamic component of displacement 


Big-end extremity of contact pattern measured 
parallel to roll axis from the midpoint of the 
effective length 


Small and extremity of contact pattern measured 
parallel to roller axis from the midpoint of the 
effective length 


Maximum permissible distance of big-end pattern 
extremity from midpoint of effective length 
measured along race 

Maximum permissible distance of small-end pattern 
extremity from midpoint of effective length 
measured along race 


Impedance matrix = Ky +1 VBy 


Other notations as defined in text. 


in. 


in. 


in. 


GREEK SYMBOLS 
Angle between roller axis and line of action of 
flange reaction 
Outer ring contact angle 
qd, cosB/E, 
dq cos(8-t)/E,. 
Displacement 


Linear displacement in x direction 


Linear displacement in y direction 


Linear displacement in z direction 


Approach of inner race to outer race at midpoint 
of effective length 


Approach of roller to race at x 


Approach of roller to outer race (cup) at qth 
roller 


Approach of roller to inner race (cone) at qth 
roller 


Residues of simultaneous equations 


oe v2) 
Roller elastic constant = ee 
E 
401 - ve) 
Race elastic constant = —————— 
ER 


Angular rotation about x axis 
Angular rotation about y axis 
Angular rotation about z axis 
Frequency cz rotation 


Poisson's ratio for roller 


xi 
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radians 


radians 


fe) 
radians, 

° 
radians, 

, ° 
radians, 


rad/sec. 


Poisson's ratio for race 


Material density 

Included roller cone angle 

Circumferential roller position 

Angular velocity of roller about its own center 
Orbital velocity of roller 


Crown drop 


xii 


lbs/in? 
radians,° 
radians,° 
rad/sec. 
rad/sec. 


in. 


SUBSCRIPTS 

Symbol Description 

b Refers to bearing 

cg Refers to center of gravity 

E Refers to roller 

F Refers to flat 

Z Refers to gyroscopic 

i Index, i = 2, 30or i=x, y, z 

i,j Refers to index of stiffness matrix; i.e., 


force in i direction due to displacement 
in j direction 


P Refers to pedestal 

q Refers to roller circumferential position 

R Refers to roller 

T Refers to total 

x Refers to x direction 
y Refers to y direction 
z Refers to z direction | 
1 Refers to outer race 


2 : Refers to inner race 


SECTION I 


INTRODUCTION 


The original Rotor-Bearing Dynamics Design Technology Series AFAPL-TR- 
65-45 (Parts I through X) included a volume, Part IV(1), which presented 
design data for typical deep-groove and angular contact ball bearings. The 
data was presented in graphical form and consisted of direct radial stiff- 
ness, load carrying capacity, and load levels. In addition design guide- 
lines and limitat »ns were discussed. The major deficiencies of this 
original volume were that centrifugal effects due to high speed were 


ignored, and axial and angular stiffness information were omitted. 


Subsequent to the publication of Part IV, several extensive treatments of 
rolling element bearings including elastohydrodynamic, thermal, and cage 
effects have been published. The computer program of Mauriello, LaGasse, 
and Jones (2) considers both eiastohydrodynamic and cage effects for ball 
bearings. The more recent computer based design guide prepared by 
Crecelius and Pirvics (3) treats elastohydrodynamic, thermal, and cage 


effects for a system of ball and roller bearings. 


Thus, very sophisticated analytical tools are available for the design and 
application of rolling element bearings. Neither of these tools, however, 
provide the user with the stiffness matrix required for solution of rotor 
dynamics problems. In addition both computer programs are very large and 


require an extensive computer facility for use. 


Part II(4) of the revised series provided an update of the original 

Part IV(1). Those aspects of the original Part IV(1) which treated general 
design aspects of ball bearings, load capacity, speed limitations, etc. 
were deleted since their coverage is superficial compared to the more 
sophisticated computer tools now available (2,3). Only those parts 


directly connected with preparation of input for the rotordynamic response 


programs (Part 1(5) of the revised series) were retained. The stiffness 


data included in the original Part IV were also updated. 


The present volume (Part III of the revised series) extends the treatment 
of rolling element bearings to the tapered roller bearing. The complete 
stiffness matrix is calculated including centrifugal effects. Considera- 
tions such as elastohydrodynamic and cage effects are not included since 
they have little influence on the calculation of tapered roller bearing 
stiffness. The resulting program (Appendix) is reasonably small and easy 


to use. 


SECTION II 


ANALYSIS 


2.1 General Bearing Model and Coordinate System 


Accurate calculation of the lateral dynamic response of a high-speed rotor 
depends on realistic characterization of the support bearings. In the most 
general case, both linear and angular motions are restrained by the support 
bearings at the attachment location. In the analytical model, the reaction 
force and the reaction moment of each bearing are felt by the rotor through 
a single station of the rotor axis As schematically illustrated in Figure 
la, a coil spring restraining the lateral displacement and a torsion spring 
which tends to oppose an inclination are attached to the same point of the 
rotor axis. A complete description of the characteristics of the support 
bearings, however, involves much more than the specification of the two 


spring constants. This is because: 


. The lateral motion of the rotor axis is concerned with two 


displacement components and two inclination components. 


. The restraining characteristics may include cross coupling 


among various displacement/inclination coordinates. 


. The restraining force/moment may not be temporally in phase 


with the displacement/inclination. 


. The restraining characteristics of the bearing may be 
dependent on either the rotor speed or the frequency of 


vibration, or both. 


Bearing pedestal compliance may not be negligible. 
To accommodate the above considerations, the support bearing characteristics 
are described in Reference 5 by a four-degrees-of-freedom impedance matrix as 


defined in Equation (1): 


be Se art A, SCT I 


(FIG la) Bearing Stiffness Model 
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(FIG tb) Bearing Location Coordinate 
System 


Rot 7 by: Hy ) 


where Wy is a column vector containing elements which are the two lateral 
displacements (6,5 6) and the two lateral inclinations (9, 0) of the 


rotor axis at the bearing station N. 


Employing a right-handed Cartesian representation in a lateral plane as 
depicted in Figure lb, the z-axis is coincident with the spin vector of the 
rotor. The x-axis is oriented in the direction of the external static load, 
and the y-axis is perpendicular to both z and x axes forming the right-handed 
triad (x, y, 2). (6,5 6) are respectively lateral lineal displacement 
components of the rotor axis along the (x, y) directions. (0,, 0.) are 
lateral inclination components respectively in the (z-x, z-y) planes. Note 
that oy is 4a rotation about the y-axis, while 0, is a rotation about the 


negative x-axis. 


Zy is a complex (4 x 4 matrix), and in accordance with the common notation 


for stiffness and damping coefficients, may be expressed as 


zy * Ky + ivB, (2) ! 


where Ls is the stiffness matrix and B, is the damping matrix. v is the 


N 
frequency of vibration. Most commonly, lateral lineal and angular displace- 
ments do not interact with each other so that the non-vanishing portions of 


Ky and By are separate 2 x 2 matrices. That is 


lineal 


(Ky) 0 0 
0 


angular 


EW oes i 
0 . F 
Be | oo 0 (4) 
2 9 eee 


Accordingly, a total characterization of a support bearing would include 


sixteen coefficients which make up the 4 (2 x 2) matrices: : 


XX xy ' 
(9) neal ” (5) ae 


~ * 
YY | lineal 


B B : 
XX xy £ 
(B) 1 ineal (6) | = 
B B | 
y YY | lineal 
: 
K K 
XxX xy 
(D angular = (7) . 
K K 
ai YY J angular 
B B 
XX xy 
(3) angular - (8) 
B ms B 
7 yy angular 


In the event that the pedestal compliance is significant, then the effective 


support impedance can be calculated from 


yi ag) (9) 


where subscripts "p" and "b" refer to the pedestal and bearing respectively. 


Note that both pedestal inertia and damping may be included in Z,° 


2.2 General Bearing Support Characteristics 


The function of a bearing is to restrict the rotor axis to a nominal axis 
under realistic static and dynamic load environments. Deviation of any 
particular point of the rotor axis from the nominal line can be character- 
ized by three lineal and two angular displacements. These may he designa-~ 
ted as (5 oe 6 Oy» 0,) in accordance with a right-handed Cartesian 
reference system. The z-coordinate is coincident with the reference axis 
and is directed toward the spin vector. (0,5 9, are rotor axis inclina- 
tions respectively in the z-x and z-y planes. The x-coordinate is directed 
toward the predominant static load; e.g., earth gravity. Ideally, the 
bearing would resist the occurrence of any displacement so that the 
reaction force system imparted by the bearing to the rotor is generally 


expressed in matrix notation as 
| ee ee (10) 


F is a column vector comprising the five reaction components (FY Fy F 
Mis Mi) while x is the displacement vector (85 by» 5 0» 0,)- Z is 
a (5 x 5) matrix containing the elements 25, with both indices (i, j) 


ranging from 1 to 5. The values of 244 characterize how rotor displacements 


are being resisted by the bearing. 


From the standpoint of dynamic perturbation, distinction is made between 
a static equilibrium component and a dynamic perturbation component for 
both the displacements and the reactions. Thus, 

28 oe ee Oe Oe OE aE. (11) 
(x', F') are respectively presumed to be infinitesimal in comparison with 
(x0» FO: Accordingly, Ry are regarded as dependent on x but not on x’. 
To illustrate the idea of perturbation linearization, one may examine the 


one-dimensional load-displacement curve shown in Figure 2. 


ee EO IE gil a” seme re 


Figure 2. 


Linearization of Tapered Roller Bearing Stiffness 


a, 


As illustrated, the load-displacement relationship is a 10/9 power law in 
accordance with the Hertzian point contact formula. It is mot possible 
to describe the entire range by a linear approximation. However, if a 
small dynamic perturbation is taken around a static equilibrium point, 
sy < é » the small segment of the load-displacement curve can be 
approximaf ed by a local tangent line. The corresponding force increment 
is 


er » 202) 


where é'Y is the incremental displacement. OF / 36, will depend on the 


amplitude of 6 . 
[e) 


The question of history dependence is resolved by regarding x' as periodic 
motions at any frequency v of interest, and he accordingly would have both 


real and imaginary parts and may also be dependent on both the rotor speed w 


and the vibration frequency v. 


To avoid notational clumsiness, the primes will be dropped from (F’, x") 
which are understood to be dynamic perturbation quantities unless the sub- 


script "o" is used to designate the static equilibrium condition. 


2.3 Tapered Roller Bearing Characterization 


In many ways the tapered roller bearing is much simpler to model from a rotor 
dynamic point of view than a fluid film bearing. In gencral, the following 


two simplifications can be made: 


. The restraining characteristics do not include cross coupling 


among the various displacement/inclination coordinates. 


The restraining force/moment is normally temporally in phase 


with the displacement/inclination. 


& 


Figure 3 shows a tapered roller bearing referred to in an orthogonal xyz ' 


Jee, 


coordinate system. The outer ring is fixed but the inner ring may move 
with respect to the coordinate system. Both rings are free to rotate about 


their axes. 


Three lineal displacements, bys oo So and two angular displacements, oN 

Oy are required to define the spatial position and attitude of the inner | 
ring when it is displaced from its initial position. For purposes of deriva- 
tion the initial situation is that existing when the bearing's end play is 
just taken up in the thrust direction. Figure 4 shows these displacements 

in the positive sense. Figure 5 establishes the convention of the roller- 


position index q. 


2.3.1 Stiffness 


The total characterization of a tapered roller bearing's stiffness 


can be expressed by the matrix. 


lor oF aF 
—s —* — + 

3x ay az 

oF oF oF 
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Figure 5. Tapered Roller Bearing Index, q 
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The lineal and angular stiffness matrices (Equations 5 and 7) can 


be derived from Equation (13). For example: 


oF oF 
—* eee. 
ox ay 
(14) 
(OT ineat j 
oF aF ; 
— a: 
ax ay 
aM aM | 
wees —_x 
30 30 
x y 
(K) = (15) 
=angular 
oM aM 
aD & 
30 a0 
x y 


Note that although the axial components of stiffness are not utilized 

by the lateral rotor dynamics program (5), they have been retained in 

the general tapered roller bearing stiffness matrix, Equation (13). ‘he 
axial stiffness would be required, for example, if the reader was calculat- 


ing the axial natural frequency of a tapered roller bearing mounted shaft. 


2.3.2 Damping 
An extensive search of the literature revealed no experimental damping 
data for tapered roller bearings. As the current state-of-the-art 


does not permit accurate calculation of tapered roller bearing damping, 


no damping data is included in this report. 


2.4 Tapered Roller Bearing Under Combined Loading 


Solution for the stiffness matrix of a tapered roller bearing under com- 
bined loading is a tedious problem and requires the use of a digital 
computer. In this section, the derivation of the solution is described. 


A computer program for obtaining the solution is included in the Appendix. 


° 2.4.1 Bearing Applied Forces and Moments 


As the result of the five displacements described previously in 
. Figures 3 and 4, there are the reactions Fye Fy Fi» and My and 
J t J ' ' 
My Fe Fy and Fo are forces. My and My are moments. All are 
shown in their positive sense in Figure 4, External forces Fy and 
F. may be applied at the inner ring center. The senses of the 


' 
signs are the same as for the reactions Fe and Fie 


2.4.2 Koller Geometry 


t 
Figure 6 shows the boundary dimensions of a typical tapered roller. 


Roller mass, moment of inertia, and location of the center of 
gravity are calculated assuming the roller is a flat-ended, trun- 


cated cone bounded by R Ro» and hee 


1° 
In general, the big-end face of the roller is not flat but is a 
sphere having the radius R, which is generally a proportion of the 
slant height, hes of the untruncated roller cone. Roller crown and 
corner breaks are also omitted from mass and moment of inertia 


calculations as their contributions are second order. 


Figure 7 is a more complete sketch of the roller showing the details 
of the roller crown. The big-end spherical surface is neglected 


here also. 


t is the included angle of the roller cone and is obtained by itera- 


tion of 


-1 ,d tT 
tan 3 sin (8B - > (16) 
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Figure 7. 
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From Figure 7 
2 F,\2 


where R. is the crown radius and he the length of the flat portion 
of the roller profile. In a fully crowned roller, the flat length 


is zero. 


Qe is the effective length of the roller load-carrying surface. 
The actual working length for any loading must lie within Loe B 
and By are the corner breaks. Their shapes are unimportant as long 


as they blend smoothly into the crowned surface. 


Vv is the drop of the crown and is measured at the extremes of the 


effective length of the roller. 


L 
v= H-YR2- CS? (18) 
Hy = —t (19) 
2tan(>) 
he T T 
Hy = Ho - Zz cos >) + Vsin(>) ~ By (20) 
- T T 
Hy = Hy + > cos >) + Vsin(>) + Bo (21) 
T 
Ry = H, tan(>) (22) 
R, = H,tan(S) 23 
2 * gtan 7 (23) 


Let the cone corresponding to Hy have the mass my) and a moment of 
inertia about its center of gravity aT . 
cg 

Let the cone corresponding to Hy have the mass 0) and a moment of 


inertia about its center of gravity I 


2° 
ce 


where p is the material density in 1b/in?, 


Then the distance X from the big end of the roller at H, to the 
center of gravity of the roller is X' 
mH 3H 
2:2 1 
(Hy - >) 
ys (28) 


The moment of inertia Lee of the tapered roller about its center of 
gravity at X is 


H 3H 
“ a2 Ry ~— _ Ry? 
a els Tar ae 5 a z Xx") 


(29) 


Later the distance x, being the distance left from Hy to the center 
of gravity of the roller, will be required. 


27 Hy - x’ (30) 


The slant height, hee of the truncated roller cone is shown in 
Figure 6 and is 


2 
sir G 


g = 


(31) 


and the big-end spherical radius, Re» is a proportion of Lee 
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V is the radius from the roller centerline to the point of contact 
of the roller and inner-race guide flange and the flange reaction, 


It is directed at an angle, a, where 
- Vv 
a ® sin ~ (=) (32) 


The lever arm of the flange reaction about the midpoint of the work- 


ing surface of the roller at Hy is 


2 = | RE - Ry - Hy + | sina ; (33) 


Figure 8 is an enlarged view of the race profile showing the crown 
drop VY which is measured at a distance G from the end of the 
effective length. The contour is the same at both ends of the roll. 


If the radius Ro is known, the drop at G is 


e . 2 <_< 2 
yrs a -\R-G G) (34) 


If the drop is known and the radius R. is not, the radius is 


(Z-6?-chr-w7) 2 


e 2 
c Sor Ber St 


(35) 


2.4.3 Roller Equilibrium 
Figure 9 shows the forces and moments acting on a roller which is 
in contact with both outer and inner races and with the inner ring 


guide flange. 


In the following discussion, the subscripts 1 and 2 refer to the 


outer and inner contacts, respectively. 


Py and Po are the contact loads. M) and My are contact moments 
resulting from nonuniform loading along the roller'’s length. Fo 


is the centrifugal force and Mo is the gyroscopic moment. The 
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latter acts at the center of gravity of the roller which is located 


the distance X from the central plane of the roller which contains an) 


the midpoint of the effective length. 


The centrifugal force and the gyroscopic moment are 


Es T 4,2 
F. = (m, + mo) {5 + X sin(8 - 2 he (36) 
. 
M. TpMpupsin(s = 35) (37) 


where X, is the orbital velocity of the roller and we the angular 


velocity of the roller about its own center, both in radians/second. 


1 dcos(8 - 5) dcos(8 = 3) 
% 8 F At od + 24 F - | (38) 


dcos(B - 5) 2 
aa = 


E 


Qy and 25 are the input angular velocities of outer and inner rings 
in radians/second. P, is the reaction of the inner-ring flange on 


the roller. 


In the present problem, we are concerned with external forces applied 
to the bearing inner ring along x and/or z (Figure 4) only. There 
may also be initial linear displacements along any or all of the 
coordinate axes, x, y, and z; and initial rotations about x and y. 
These initial displacenents do not change when external forces are 
applied along x and/or z. However, when initial rotations are pre- 
sent about x or z, cperating displacements may occur along x and/or 

y as the case may be. The system, therefore, has the possibility 

of three degrees of freedom; i.e., working linear displacements along 
any or all of the axes x, y, and z. If initial displacements exist 


about x or y, working displacements in these modes are prevented. 


23 


The approach of the inner race to the outer race along the line 
defined by 8 for a roller at azimuth ¢ is 


P 
D 
= Ww ie wt wy. 7 
A (6, + Ss) sinb + {(6, + 84) cosd + (S, + 6 sing Fs! 


cosB + S{Esing + dsin(Z)}{(o, + o")sing + (0, + 0")cosg} 
(40) 


Py is the diametral clearance or the total diametral play of the 


inner ring relative to the outer ring before loading. 


The azimuth angle, 9, is related to the roller position index, q, 
through 


o = cb (41) 


n 


where n is the number of rollers. 


The double-primed items in Equation (40) are the initial displace- 


ments in the several modes. 


Also, as a result of the initial misalignments which may exist about 
x and/or y, the inner race at the qth roller may be misaligned the 


amount @Q. 
Qe = (o, + oY) sing + (9, + 9) cos¢ (42) 


If Ay is the approach of the roller to the midpoint of the outer 


race, the approach Ay of the inner ring to the roller at its midpoint 


is 


(A - b,)cos(a - 4) 
A = ee nce (43) 


4 cos(a + 3) 


If 0, is the misalignment of the roller relative to the outer race, 


the misalignment 85 of the inner race relative to the roller is 


0, = 0-% (44) 


Misalignment is positive if it tends to squeeze the big end of the 
roller more than the little end when the big end is at the left. 


Figure 10 illustrates the geometric intersection of a roller and 


raceway. 


The profiles of race and roller bodies are referred to an XY 
coordinate system. Note that the X axis is positive to the left 
of the origin. 


The equation cf the race surface is 


Y= 0 (45) 


The equation of the flat portion of the roller or the element of the 
basic roller cone is 


Y = Ay + Xrano, (46) 


The equation of the crowned portion of the roller profile is 


2 2 2 
(X - Hsino,) + (Y + Heos0, - ) Ro (47) 


The subscript i is 1 for an outer contact and 2 for an inner contact. 


The intersections of the race and the crowned roller surface occur at 


xX and X 
i Bs 


{ 2 2 
“a Ro - (Heosa,; - 4,) + HsinO, (48) 
~\fr2 - (Hcos9, - 4;) + Hsin, (49) 


% 


i 
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Xa; and Xp, must be within the projected extremities of the roller 
crown. That is 


* 
X, < *a, (50) 
m 2M 6» 
oo RE 
where 
k he 
X, . > cosG, + Vsino, (52) 
i 
x 4. 
x 2 z 6080, + VsinO, (53) 
i 


If the quantity under the radical in Equations (48) and (49) is zero 


or negative, there is no contact between roller and race. 


g 
If a cos0, > X, , there is also no contact. 
2 i- A, 
If X, > x" Xe ed t lt x" 
; s set equal to . 
a Ay A; Ay 


If X% < a ’ 6 is set equal to % 


i i 
he Le te 
lf > cos, > ‘3, >- > cos, and ‘s > os cos, : 
as 
the value of X, is X, Barry ae (54) 
i i 5 
From Figure 9 the conditions for roller force equilibrium are 
Tt 
“Py cosB + P, cos(6-t) - P, sin(B - a 41) + Fo 0) (55) 
Tt 
Py sinB + Py sin(p-t) + P, cos(B ~ ig a) 0 (56) 


Equations (55) and (56) are a set of simultaneous nonlinear equations 


in which the variables are 4) and 8) at the outer contact of the 
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particular roller. 


The flange reaction P, is obtained by taking moments about the roller 


midpoint. 
(-M, +M, - M+ PF. Xeos(p - 4) - 2 (Pp, - P,)sin()} 
Pp = 1 2 G c 2 2 1 2 2 
3 g 
(57) 
From Figure 10 the intrusion of the roller into the race is 
hp 
S, = A, + Xtano, Ix| < > coso, (58) 
. he 
A, = \R - (X-Hsino, ) - Hcos0, + A, [x| > > cos0,; 
(59) 


The derivatives of AL with respect to 0, will be required later and 


are 
dA £ 
ee ee |x| < = cos8, (60) 
d& 2 - 2 i 
i cos 0 
i 
(X-HsinO, )Hcosd L 
OO ig. Re A Deh |x| > = cosO, (61) 


st eee L 
do |] 2 2 
i Ro - (X-HsinO,) 


Lundberg (6) gives the approach Ay of two cylindrical bodies pressed 


together with the uniform loading P, as 


(ng + nQ) Ki =X, 
R E A 
A. sa me, Toons {1.8864 + &n (2p (62) 
Ne and Np are elastic constants for race and roller, respectively, 


having the form 


= 4Q - v?) (63) 


Rid 
R,E ERE 


where v is Poisson's Ratio and E is the modulus of elasticity. 


by is the semi-width of the pressure area in the rolling direction, 


(np + ne) 1/2 
by e [a= Py qd, q+ Si YD) (64) 


C; is 1 for 1 = 1, corresponding to an outer contact; and ~1 for 


i = 2, corresponding to an inner contact, 


d .cosB 
x 
d cos (B-t) 
ae 666) 
x 
where z 
d+ 2Xsin(>) 
(eee (67) 
* s(3) 
cos (5 
E. = E + 2Xsin8 + dcos(6 - ) - dcosB (68) 


The value of Py corresponding to Ay is required. This cannot be 
obtained from Equation (62) in closed form. It can be obtained 


numerically in the following manner. 


Let Py. be an estimate of Py: A good starting value is 


5 x 10° aa? 


Pye 1/9 (69) 
(X,~X3) 


An improved value of Py is 


(Ay - we) 
aia (70) 
x x dA} /dp, 


Ay is the approach of race and roller bodies calculated for the 


current estimate of Px using Equation (62). 
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da’ /dpy is obtained from Equations (62) and (64) using the current 


estimate Pk and is 


da‘ (ng + n,Q) Xx, -X 
x R E A B 
dp? on {1.3864 + gn ¢ re 


)} (71) 


Iteration of Equation (70) yields Pp, to any desired accuracy. 


The contact force, P, and the moment, M, are 


Xx 
Ad 
Po = j Pi dXx (72) 
Xp 
i 
Xa. 
1 
M, =f X P,dX (73) 
By 


Equations (55) and (56) may now be solved for Ay and 04° the dis- 
placements at the outer contact. Again, a closed-form solution 


cannot be obtained and numerical techniques are employed. 


If estimates are made of the variables Ay and Oy Equations (55) 
and (56) may not be satisfied and there will be the residues ey 
and Eo for Equations (55) and (56), respectively. Differentiating 
Equations (55) and (56) gives: 


de, dP, dP, dA Z dP, 
= -cosB ~— + cos($-t) —— sin(B - ~- a) == 
da) dA, ds. dd, 2 day 
(74) 
de dP dP, dO dP 
1 1 2 Tt 3 
ax = -cos8 ~— + cos(B-t) == 7 — Sa) 
d9, do, do, do 2 do, 
(75) 
de, dP, ; dP, dA ee dP, 
a Sine a tb Sint ht) a Se C08 (8 SS 
dé, da, dé, ds 5] da 5 


de, dP dP, do 


= sing = + pint poe ae 
Oy 


do, do, do, 


From Equations (43) and (44) 


T 
dA, . -cos(a - >) 
dh) cos(a + =) 
i er 
do, 


And, from Equation (57) 


ee : dA, da, dA, 2 dA, 
ds, 2 
dP, . de, 05 01 2 do, 
dO) 2 


+ cos(B - z- a) 


dP, dA 
2.42 T 
~- —* —)sin(D 
dA, dA, 2 


dP, 


do, 


If Ay and oH are current estimates, improved estimates are: 


: de, 
1 do, 
de, 
2 do 
A, = At - : 
1 1 de, de, 
dA, de, 
de, de 
da, dO, 


31 


(77) 


(78) 


(79) 


(80) 


(81) 


(82) 


es 
dA, 1 
be : 
di, “2 ; 
@, = Oo! - (83) 
1 1 de, ao 
dA, do, 
dA, do, 


The determinants in Equations (82) and (83) are calculated at current 


estimates. 


The derivatives of Py and M; with respect to A; and 05 are | 


aP, Ka, dp, dd 
dA -[ aa. aa, %* Pe 
i B x i 
i ‘ 
{ 
| 
dP X dp. dA | 
ane ee Ay a —* ax (85) 
do dA. do, 
i X3 x 
di 
aM, Xa dp. dd 
i = 4 x x : 
dA, ii Maa ds p ASE) 
“B x 
i 
dM. Ce dp. da 
1 1 x x 
— = X = = dX (87) 
d0,. “a da. do. 


The value of dp, /dd, is obtained from Equation (71) and the value of 
da /dd, is unity. 


If Equations (43), (44), (55), and (56) are differentiated with 


respect to A, there results four equations which are linear in 


dA, /dd, dd,/dd, de, /da, and d0,/dd and from which all four deriva- 
tives can be obtained. Of the four derivatives, only dA, /ds and 


do, /da are of interest here. 


dP, - dP ds, dP, 
[-coss a, - sin(p - 77 a) a, aA + [cos¢e-t) aa, - 


¢ dP, da, dP) a dP 
sin(S - $-a) apo) ay + [-cos8 ag - sine - F- ©) ae] 
2 1 1 
do dP dP, , d9 
1 2 t 3 2 
— + [cos(e-z) ®, - sin(B - 7 - a) a. ae = 0 (88) 


dP ? dP, da, dP. 
[-sins aa. + cos(B - >7 a) =| in + [ since-z) ar + 
1 1 2 
: dP, dA, a, : 2] 
cos(8 - z7 a) i, in + [-sine a, + cos(B - 77 a) do, 
do dP dP d9 
1 2 T 3 2 = 
aia + [ since-r) a0, + cos(f - ies a) d0, a 0 (89) 
( aS 
bra rare a ae 1 (90) 
da t, dA 
cos(a - >) 
da da 
: 1 2 
a6 ae (91) 


Equations (88) through (91) are easily solved for dA, /da and de /da. 
dA ,/do and do, /d9o are obtained in a similar manner. 


2.4.4 Bearing Equilibrium 


The reactions of the bearing on the shaft at the central plane of the 


roller are 
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n - 
Fis P 92 
‘ cose L, 1,284, (92) 
q 
n 
F' = cos P;} sin (43) 
y cos x, 1, +4 
n 
Fi = sing Py (94) 
q=l 4 
s 1 Tt 
Mos Sy [ F(esins + dsin(Z)}P) + My ]sine (95) 
* q=1 ~~ q q q 
te 22: i PaaS 
My = &, [ desing + dein) Pr. + my Js084, (96) 


Considering the three-degree-of-freedom system, the inner ring is 
acted upon by the external forces FS and Ee and may have working 


displacements along x, y, and z. Equilibrium requires that 


Fe + EF, = 0 (97) 
Fy = 0 (98) 
Fo + F = 0 (99) : 


Here the variables are bye Sys and 5. Again, a direct solution is 
not possible and numerical methods must be employed. 

For initial estimates bye by 5 of the variables Equations (97), 
(98), and (99) may not be satisfied and there remain the residues 


ey E4, and 3° Improved values of the variables are 


The right members 


current estimates 
d(é,, bys 6) 
d(4., bys 6) 


d(4,, bys 6) 


(104) 


(105) 


(106) 


Although only the above derivatives are required in determining the 


equilibrium of the system, the complete matrix is required for 


stiffness calculations. 


dF’ n ori. 
Ca aC) = cosB x, cost, T5808) (107) 
dF n ae 
Wit) = cosB x cosy, FEB 05 0,) (108) 
qF? n “ae 
F559 89839 9 Oy) Teen x 45696710, 9,) see 
a = Liaegeued fies ae 
TE 558, = S = Li sinBt+d sin(>) T5558, 0) * 
“, 
TCR CRD sing, (110) 
aM! nr) . ma 
T5580) = ple sin8+d sin(5) HE 58,80) + 
amy | 
Ce ] cose, (111) 
where 
dP, dP, dt, dP) do, an 
ae: ere ne Pate, || aren arte, comma 
dB 8518, 9,50,) la dq * do, | FY 8598 50,505) 


(6,6 ,6.,0,0 -| 77 
x? y’ z’ x? ) dA, aN 


The derivatives of A 


dM, dM) dA 


q 
= dé, 
ds, dO 
q 


cosp cos 
%q 
cosB si 
s ne, 
sing 
Lig sinf+d sin(=)}sing 
2 2 q 


1 T 
={E sint+ ~)} 
DIE sin3+d sin(>) cos}, 


2.4. Effect of Unloaded Roller 
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+ 


dA 


ar 
re ae AE Hed 


y’ 2’ x? y 


d@ 


fe a ae Sa 
Co a 


(113) 


and o with respect to the inner-ring displace- 
ments are, from Equations (40) and (42) 


(114) 


(115) 


(116) 


(117) 


(118) 


(119) 


(120) 


In some instances one or more rollers may be out of contact with the 
inner race while in contact with the outer race and the inner-ring 


The conditions for equilibrium of such rollers are 


~P, cosB - P,sin(8 - 3 - a) + F 0 (121) 
~P) sing + P,cos(8 - >: a) = 0 (122) 
where 
“M,-M, - 3 P,dsin(>) + F, X cos(B - 3) 


g 


Here the variables are Ay and 0). Initial estimates Ay and oF will 
generally fail to satisfy Equations (121) and (122), and there will 


be the residues ey and Eo: 


Improved values are 


, de, 
1 dO, 
de. 
ae, do, 
A = A! — (124) 
1 dey de, 
da, do, 
de, de, 
dA dO, 
de 
irene € 
da, 1 
de, : 
dA, 2 
) = Q! - (125) 
1 1 de, de, 
dA, doy 
de. de 
( 
day dO 


The right members of Equations (124) and (125) are evaluated at 
current estimates. Iteration of Equations (124) and (125) yield 


ay and OF to any desired accuracy. 


The derivatives required in Equations (124) and (125) are 


where 


Rollers which 


in evaluating 


nothing to the stiffness matrix since P 


dP, e dP, 
= cosg an sin(g - 77 a) a” (126) 
1 1 
dP dP 
a sg Ne , 
- cosg q6 sin(p 7) a) FTG) (127) 
1 1 
dP dP 
~ sing eee cos(t - + - a) = (127) 
dA 2 da 
1 1 
dP - dP. 
- sinB qo. + cos(B - 77 a) qo, (128) 
dM P 
1 1 T 1 
~ 72 dsin(3) ds, 
(129) 
R 
dM dP 
1 1 T 1 
- a0, a5 dsin(>) q, 
ee (130) 
g 


are out of contact with the inner race must be considered 
the bearing's reactions. They, however, contribute 


1 and My for these rollers do 


not change with changes in the inner-ring displacements. 
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SECTION III 


APPLICATION OF COMPUTER PROGRAM 


The analysis of Section II has been programmed in Fortran IV for a 
digital computer and is suitable for use on the CDC 6600. A program 


listiug is presented in the Appendix. 


3.1 Sample Test Case 


To illustrate a typical case consider the bearing in Figure 11. This 
is a tapered roller bearing assembly modified for high speed operation. 


The geometry of this sample bearing is summarized below. 


Number of rollers 37 
Roller diameter at midpoint 2913 in, 
Pitch diameter 5.0 in. 
Contact angle at outer race 14°40" 
Effective length of roller -6001 in, 
Roller big-end spherical radius 0.8 in. 
Radius from roller centerline to point of big~ 0.75 in. 
end spherical surface with inner race flange 
Roller crown radius 100 in. 
Roller small-end corner break 02 in, 
Roller big-end corner break 03 in. 
Crown drop gage point 03 in, 


The operating conditions for the sample case are: 


Rotational speed = 20,000 rpm 
Load Condition #1 
Thrust Load 
Load Condition #2 
Thrust Load 
Radial Load 


3,000 lbs. 


3,000 lbs. 
700 lbs. 
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Figurell. Sample Tapered Roller Searing Assembly 
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3.2__Input_ Format 
Figure 12 presents the input data format and Figure 13 shows the actual 
input data for Load Conditions #1 and #2 of the sample case. 


3.3 Output Format 

Figure 14 presents the output data for Load Cendition #1. The input data 
are summarized in Figure 14, followed by the output data including the 
internal load distribution as well as various other stress and displace- 
ment parameters. The stiffness matrix is given on the last page of Figure 
14, 


The output data for Load Condition #2 are presented in Figure 15. 
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